The Role of Collective Neutrino Flavor Oscillations in Core-Collapse Supernova Shock Revival 
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We explore the effects of collective neutrino flavor oscillations due to neutrino-neutrino interactions 
on the neutrino heating behind a stalled core-collapse supernova shock. We carry out axisymmetric 
(2D) radiation-hydrodynamic core-collapse supernova simulations, tracking the first 400 ms of the 
post-core-bounce evolution in 11.2-Af0 and 15-Af0 progenitor stars. Using inputs from these 2D 
simulations, we perform neutrino flavor oscillation calculations in multi-energy single-angle and 
multi-angle single-energy approximations. Our results show that flavor conversions do not set in until 
close to or outside the stalled shock, enhancing heating by not more than a few percent in the most 
optimistic case. Consequently, we conclude that the postbounce pre-explosion dynamics of standard 
core-collapse supernovae remains unaffected by neutrino oscillations. Multi-angle effects in regions 
of high electron density can further inhibit collective oscillations, strengthening our conclusion. 
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I. INTRODUCTION 

Stars with a zero-age main-sequence (ZAMS) mass 
greater than ~(8— 10)-Mo undergo core cohapse at the 
end of their lives. When the inner core reaches nuclear 
density, the stiffening of the nuclear equation of state 
(EOS) induces core bounce, launching a strong shock 
wave that slams into the still infalling outer core. Created 
at a mass coordinate of ~O.5-M0 and endowed with the 
kinetic energy from infall, the shock initially moves out 
rapidly in mass and radius, but dissociation of heavy 
infalling nuclei and neutrino losses from neutronization 
and thermal processes in the hot region behind the shock 
sap its energy. As a result, the shock soon succumbs to 
the ram pressure of the outer core and stalls at a radius 
of -(100- 200) km. 

Finding the mechanism that robustly revives the 
stalled shock to blow up a massive star in a core-collapse 
supernova has been the primary objective of core-collapse 
supernova theory for decades, but so far success has been 
limited. The neutrino mechanism [l|,[2[, based on the net 
deposition of neutrino energy by chargcd-current absorp- 
tion in the semi-transparent gain region below the shock, 
has been shown by radiation-hydrodynamics simulations 
to work in its purest, spherically-symmetric (ID) form 
only in the lowest-mass massive stars with oxygen-neon 
cores and the most accurate ID simulations fail to 
explode more massive stars @, In axisymmetry (2D), 
neutrino-driven convection and the standing-accretion- 
shock instability (SASI, e.g., [sl-fioj) increase the efficacy 
of the neutrino mechanism (by enhancing heating |lll4l3j 
or reducing cooling by neutrinos [l3|). Neutrino-driven 
exp losions in 2D simulations have been reported [l3l [l5l . 
[la |. though only in models using the softest variant of 
the Lattimer-Swesty EOS which is disfavored by the 
recent discovery of a 2-Mq neutron star [31 : but leads to 
a compact protoneutron star (PNS) and a hard neutrino 
spectrum [ll, favorable for heating due to the 
dependence of the neutrino absorption cross section. 



Other multi-dimensional phenomena have led to pro- 
posals of alternatives to the neutrino mechanism: Rapid 
rotation in combination with rotational magnetic field 
amplification can lead to magnetorotational explosions 
with bipolar morphology [20|-|23| . But rapid core rotation 
may be p resent only in a small fraction of all massive 
stars |24l . [25| . ruling out this magnetorotational mech- 
anism for the garden-variety core-collapse supernova. 
Burrows et al. [261 , [27| proposed an acoustic mechanism in 
which non-linear PNS oscillations, driven by turbulence 
and accretion downstreams, emit sound waves into the 
region behind the shock that steepen to secondary 
shocks, injecting heat and reviving the stalled shock. 
The acoustic mechanism leads to perhaps too late, too 
weak explosions and bleeding of PNS oscillation power 
to numerically unresolved daughter modes may diminish 
its relevance in nature [28j . 

Exploratory 3D simulations [1^, HO] are suggestive of 
the possibility that the additional degree of freedom over 
2D and the physical nature of turbulence in 3D could 
render the neutrino mechanism robust. But full 3D 
neutrino radiation-hydrodynamics simulations must be 
awaited before 3D can be declared the missing piece in 
the core-collapse supernova puzzle. 

The marginality of the various proposed mechanisms 
combined with nature's robust ability to produce explo- 
sions in massive stars up to at least — 20-Mq [sTI. |32|| 
makes one wonder: Is there important physics missing 
from core-collapse supernova models? With the current 
standard set of physics included in models, we might, for 
example, miss an early hadron-quark phase transition in 
the PNS, leading to a second collapse and bounce and a 
second shock wave that revives the first. This possibility 
was proposed by [33l . [3^ , but requires a soft hadronic 
EOS that has now been ruled out. 

In this article, we consider new physics that has not 
previously been included in the core-collapse supernova 
problem: self-induced collective neutrino flavor oscilla- 
tions. Neutrino oscillations have long been known to 
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occur in vacuum {e.g., [35| l and in matter, mediated by 
neutrino-electron scattering (via the Mikheyev-Smirnov- 
Wolfenstein (MSW) effect, 0, [ll). Less appreciated, 
until recently, has been the possibility of oscillations 
occurring due to neutrino-neutrino forward scattering 
in regions of high neutrino number density. This self- 
induced oscillation process was first discussed in [s^ and 
then explored in a series of papers |39I446|. Thereafter, 
following the simulations of Duan et al. |47l. lisj , it has 
received much attention recently as a process occurring 
in the core-collapse supernova environment and leading 
to flavor conversions of neutrinos and antineutrinos of al- 
most all energies {e.g., j49l - l68j . Also see the review f69j). 

The most intriguing result of these so-called "collective 
oscillations" (flavor conversions take place collectively 
over all energies) is an almost complete exchange of 
electron neutrino {v^) and antineutrino {Vf.) spectra 
with the spectra of the heavy-lepton neutrinos and 
antineutrinos Vx G {i/f^,Vr,v^,VT}. The Vx, due to 
the absence of niuons and taus, do not interact via 
charged-current processes in core-collapse supernovae. 
They are created by thermal processes in the PNS core 
and decouple from matter at smaller radii and higher 
temperatures than and Df,. Hence, the initial Vx 
spectra are much harder. Due to the dependence of 
the charged-current absorption cross section, a swap of Vx 
and Ve/i^e spectra would dramatically enhance neutrino 
heating and may be the crucial ingredient missing in core- 
collapse supernova models, provided the swap occurs at 
sufficiently small radii in the region behind the shock 
to boost net heating. This is indeed the recent result 
obtained by Suwa et al. [tO] , where a neutrino conversion 
radii was assumed at a radius of 100 km. 

Fuller et al. [7l| were the first to propose increased 
heating due to neutrino oscillations, but their MSW 
resonance based oscillation mechanism required a large 
neutrino mass of ^ (10 — 100) eV for one of the active 
neutrinos. Akhmedov et al. (72| put forth a similar 
proposal, but both are now ruled out by stringent 
constraints on neutrino masses {e.g., [73|). 

Self-induced collective neutrino oscillations, on the 
other hand, do not require large neutrino masses, are 
a rather straightforward consequence of the Standard 
Model of particle physics, and may, in principle, occur 
in both the inverted and the normal neutrino mass 
hierarchy [gO]- In this work, we study their relevance 
for the core-collapse supernova mechanism by calculating 
approximate analytic and detailed numerical estimates 
for the radii at which collective oscillations set in 
and could influence neutrino heating. We base these 
calculations on neutrino radiation fields from 2D neutrino 
radiation-hydrodynamic simulations of the postbounce 
core-collapse supernova evolution in 11.2-Mq and 15- A/© 
progenitor stars, representative of the progenitors of 
standard Type-II supernovae. 

We find that the calculated oscillation radii, while 
reaching average shock radii, due not penetrate deeply 
into the heating region. Large shock excursions due to 



the SASI reach and surpass the radius at which oscilla- 
tions set in, but the region in which the vast majority of 
net heating occurs remains always at least ^100 km be- 
low the oscillation radius, even in the low-mass 11.2 -Mp, 
progenitor. Recent results of Chakraborty et al. [tJ, [75| . 
obtained on the basis of ID simulations, show the 
suppression of collective oscillations by very high electron 
number density in dense matter. In our 2D models, we 
find that such a suppression may be significantly weaker 
than reported by Chakraborty et al.. Nonetheless, 
the flavor conversion still happens too far out in the 
supernova, and we conclude that collective neutrino os- 
cillations do not have a significant effect on the explosion 
mechanism of core-collapse supernovae in progenitors in 
and above the explored mass range. 

The structure of this article is as follows. In Section 
Hn we review neutrino heating in core-collapse supernovae 
and in Section lllll we introduce collective neutrino oscil- 
lations and present an approximate analytic prescription 
that can be used to determine the radius at which 
neutrinos will begin to collectively oscillate. In Section 
ITVl we go on to discuss our 2D radiation-hydrodynamic 
postbounce core-collapse supernova simulations and con- 
trast the evolutions of their characteristic radii with the 
analytic estimates and detailed numerical results for the 
oscillation radius. This allows us to ascertain the impor- 
tance of collective neutrino oscillations for shock revival. 
In Section |Vl we summarize our results and conclude. 

II. NEUTRINO HEATING IN CORE-COLLAPSE 
SUPERNOVAE 

To elucidate the basics of the neutrino mechanism, we 
make a number of simplifying assumptions, which we lay 
out in the following. We assume a spherically symmetric 
mass distribution and expect neutrinos to stream freely 
outside their energy-averaged neutrinospheres. We define 
the neutrinosphere radii for each neutrino species via a 
Rosscland mean neutrino optical depth, 

and set the neutrinosphere radius R^^ = R{TjiM,i^i = 
2/3). In this expression for the energy-averaged optical 
depth, K,^. is the sum of the absorption and scattering 
opacities, and J^- is the energy density. 

The energy-dependent optical depth is a quadratic 
function of the neutrino energy, the energy-dependent 
neutrinospheres move outward with neutrino energy (cf. 
Fig. 13 of Izll)- The energy-averaged variant of the 
optical depth that we use produces a neutrinosphere that 
matches rather well the radius of the neutrinosphere of 
the average neutrino energy for each species. In the 
following, we will consider only the neutrinosphere 
radius Ri,^, since Ve decouple from matter furthest out, 
followed first by De, then by i^x. The Vx are not involved 
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in charged-current interactions, hence have the largest 
mean free path. The interact with the less abundant 
proton and, hence, also have a greater mean free path 
than the v^. This decoupling hierarchy is also present 
in the mean and mean-squared neutrino energies, giving 
{ev^) > (ep,) > (e,.,) and (e^J > {elj > (e^J, reflecting 
the fact that the least interacting neutrino decouples at 
the smallest radius and the highest matter temperature. 
Typical values for the mean neutrino energies, (cj/;), are 
~(10 - 15) MeV for j^and De and ~(15 - 20) MeV for 
(see, e.g., [H, [z3, Izll)- In the very early postbounce 
phase of a core-collapse supernova, R^^ is typically 
around ^(70 — 80) km, reaching ~(30 — 40)km within 
(200 — 300) ms of core bounce. 

We make the assumption that beyond our nominal 
neutrinosphere R^^, the radiation fields of all neutrinos 
are freely streaming with a luminosity C^. , an average 
energy of (e^.), a mean squared energy of (e^.), and 
a total number luminosity of Af^- = AttR^ , where 
$1/; is the neutrino number flux of species i at radius 
_Ry^ . We normalize the spectral neutrino distribution 
function d^^./de at the z/g neutrinosphere according to 
4:TtRI^ ded^^Jde = AttRI^^^^. The neutrino energy 
distribution function is then simply ed^^^/de. 

Outside of R^^ , conditions arise where absorption of 
i>e and via charged-current interactions with neutrons 
and protons inject more energy into the matter than is 
lost due to thermal emission and electron and positron 
capture on neutrons and protons. This is the case 
in the gain region (79j . which, in our simplified ID 
picture, extends from the gain radius Rg, where neutrino 
heating balances cooling, to the shock radius Rg. Typical 
values of Rg and Rs during the accretion phase of 
core-collapse supernovae are '--^lOOkm, and ^200 km, 
respectively [H, 111, [zS HI . 

In the gain region, the heating due to the charged- 
current neutrino-matter interactions is given by the 
absorption cross section ai,-{e^-) convolved with the 
spectral energy fluxes of both v^s and i?eS incident on 
the gain region from below, 

drAnr^nJ dea,.(e)i-^, (2) 

where is the local number density of nucleons relevant 
for interactions with neutrino species i. For Eq. we 
approximate cTuiie) = cto(1 + 3g\)/{Amlc^) x = ae^, 
where cto ~1-76 x 10"'*'' cm^ is the reference weak inter- 
action cross section, and gA is the axial-vector coupling 
constant. Note that the units of a are cm^/MeV^. This 
approximation of the neutrino absorption cross section, 
in combination with our free-streaming assumption, leads 
to a gross heating rate, 

H ~ y^(T(e^.)£^^ / drui 

~ a[{el^)Ci,^CN + {e\)Cy^cp\ , (3) 



where cn and cp are the target nucleon column densities 
(in #/cm^) for and D^, respectively. Note that Cy. 
and {el) will typically vary somewhat across the gain 
region. We neglect this variation for simplicity. Eq. ^ 
gives the integrated gross heating rate, %, but it is 
important to note that the majority of the heating occurs 
very near the gain radius due to the strong dependence 
of the rest-mass density on radius, p oc r~^ (as pointed 
out by (79I and seen in simulations. See, 6.17., Fig. 16 
of [231, and Fig. [S] of this work). Also, while neutrino 
heating dominates over cooling in the gain region, the 
latter is still significant and must not be neglected. The 
net heating rate (heating minus cooling) is estimated by 
Janka [t^ to be ?^not = 'H— C«a'H/2, which is in rough 
agreement with what one finds in simulations. 

The main take-away message from this section is that 
the heating rate % depends on the and De spectral 
fluxes as 

Our aim in this study is to explore if flavor conversions 
due to collective neutrino oscillations can increase the 
quantity >C,y^,i7^(eJ p ) to boost the heating mechanism. 
There are two aspects to this question - (i) What is the 
typical degree of flavor conversion? (ii) Does this flavor 
conversion take place at sufficiently small radii to have 
a significant effect on the total neutrino heating? To 
address these questions, we discuss collective oscillations 
in the next section. 

III. COLLECTIVE NEUTRINO FLAVOR 
CONVERSION 

A. Equations of Motion 

Neutrinos with masses mi, 7712, m^, are related to three 
flavor states v,., v^, Vr. These flavor states oscillate 
from one to another as a function of time, depending 
on the mass-square differences Am?- = m? — m? and 
mixing angles 9ij, where the indices run over (1, 2, 3). 
For an introduction to neutrino oscillation physics, see 
e.g.. Chapter 3 of [8l|. As is usual in particle physics, 
we will use units in which the speed of light c and Planck's 
constant h are equal to 1. 

In the core-collapse supernova context, a two neu- 
trino fiavor approximation is often appropriate because 
Am^i ^ |Am§i| and 6*13 ^1. In that case, the 
oscillations arc primarily between Vf. and the linear com- 
bination (v^ — Vr)l\f2, while the other linear combination 
is decoupled from the system. Our numerical results have 
been obtained with a full three-fiavor code. The 
and Ve -H- i^a; oscillations are governed by Am^^j^^ k, 
Am3i (the subscript "atm" is used, because this quantity 
determines the oscillations of neutrinos created in the 
atmosphere) and the mixing angle ^13. We take the 
absolute value of Am^^^^^ to be 2.6 x 10~^ eV^, close to 
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the best fit of experimental data [82| , while the sign is not 
known. We will use a benchmark value of 6*13 = 0.001, 
consistent with the upper limit sin^ ^13 < 0.035 [s^ . 

Neutrino oscillations can alter the flux differences 
d^i,^ /de~ d^u^/ de and /de — d^^^ /de at a given en- 
ergy. We will represent the flux differences at each energy 
e by a polarization vector P and P in a three-dimensional 
flavor space, as in [i^. At the ncutrinosphcre, all 
neutrinos are emitted as flavor states. Thus the initial 
polarization vectors are aligned with the z-direction 



P(6) 
P(6) 



d<^j^Jde - d<^i^Jde 



4$, 



(5) 
(6) 



where a vertically upward vector represents a Ve (or i>e) 
excess and a downward vector shows a excess. When 
the fluxes of i>e (or 9^) and Vx are equal, the polarization 
vector vanishes. Other directions represent states that 
are coherent superpositions of the two pure states which 
will be generated by flavor oscillations. 

In the treatment of neutrino oscillations it is important 
to consider the angle i^r^ at which neutrinos are emitted 
from the neutrinosphere at R^^ . We will therefore attach 
a label u to polarization vectors signifying the direction 
of emission, i.e., P(e, u), where u = sin^ d^j^ . Of course, 

/g duP{e, u) = P(e), and similarly for antineutrinos. 

In vacuum, the Ve (or De) oscillate to i>x and back with 
a frequency [35[ 



0.65 km" 



10 McV 



(7) 

(8) 



solely under the action of neutrino masses and mixing. 
The vacuum Hamiltonian is represented by a vector 
u;(e)B with B = ^(sin 26*13, 0, cos 2^13) where the mi- 
nus sign applies for normal neutrino mass hierarchy 
(Am^jj^ > 0) and the plus sign for inverted neutrino 
mass hierarchy (Am^^^ < 0). 

In the core-collapse supernova environment there are 
additional sources of flavor oscillation, i.e., weak inter- 



actions with the stellar material, and weak interactions 
between the neutrinos themselves. In the limit that 
neutrinos are free streaming, i.e., outside the neutri- 
nospheres, elastic scattering with electrons is the only 
neutrino-matter process that is relevant. Of these, the 
forward scattering amplitudes add coherently with the 
free propagator to introduce a potential in the flavor 
evolution Hamiltonian. This contribution due to the net 
local number density of electrons in the medium is known 
as the MSW potential lli,!!^. 



A(r) = x/2GFK(r)-ne+(r)] 



6.6 X 10^ knr 



ne{r) - ne+{r) 



1033/c 



(9) 
, (10) 



where, in the second line, we have used the value of the 
Fermi constant Gp = 1.16 x 10"^ GcY~^. The matter- 
induced contribution to the Hamiltonian is represented 
by a vector A(r)L, where L = (0, 0, 1). 

The above two contributions arc well-known, and lead 
to the traditional paradigm of core-collapse supernova 
neutrino oscillations based on vacuum oscillations and 
matter-induced oscillations [s^l- However, near the 
neutrinosphere, neutrino densities are very high, so in 
addition to ordinary neutrino oscillations due to uj{e) 
and matter oscillations due to A(r), one has appreciable 
forward scattering of neutrinos and antineutrinos off each 
other This leads to another potential, induced by 

all other neutrinos and antineutrinos, whose value at the 
neutrinosphere R,^^ is given by [3S 



= 1.1 X 10^ km 
■(10km)2 



(11) 
(12) 



E 



10^2 erg/s 



lOMeV 



The exact quantity that appears here, i.e., $^,,y = 
$1/^ -1-4$^^ , depends on our chosen normalization for 
the polarization vectors - Eqs. ([5]) and (jS]) have the same 
quantity in the denominator. Altogether, the equations 
of motion become [52| 



dP(e, u) uj{e)B x P(e, u) A(r)L x P(e, u) 



dr 



dr 



Vr{u,r) 

e ^2 



du' / de 



Vr{u,r) 

' ,P{e',u')-Pie',u') 



Vr{u',r) 



dP{e,u) _ w(e)BxP(e,w) A(r)LxP(e,u) 



Vr{u,r) 



Vr(u,r) 

Vr[u',r) 



(13) 



P(e,^) 

Vriu,r) 



P(e,») 

Vr{u, r) 



de' (P(eVP(^'))) xP(6,^) 



(14) 



^ de'{P{e')-P{e'))jxP{e, 
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Note that neutrinos emitted at angle -dn^ 
velocity, 



have a radial 



Rl 

Vr{u, r) = cosi9(u, r) ^ d 1 — u — ^ 



(15) 



and their flavor evolution has been projected onto the 
radial direction. 

To keep our discussion simple, we will often use an 
effective spherically symmetric description proposed in 
[s^ . where all neutrinos are assumed to be emitted at 
45° to the nominal neutrinosphere at R^,^ . This is often 
referred to as the single-angle approximation. We will see 
that there is no flavor change close to the neutrinosphere, 
thus this choice of a common neutrinosphere merely acts 
as a boundary condition where we specify our initial 
states. In this approximation, all neutrinos have a radial 
velocity 



mimicked by setting the mixing angle 6'i3 to a small value, 
and removing A(r) from the equations of motion. We 
use this result without proof. Now, adding Eq. ()19|) and 
Eq. (PO]) and integrating over all energies, one finds that 
the vector de(P + P) — ffZitg acts like a pendulum 
with the energy ^ 



E 



B 



de (P 



(21) 



where uj is the the average of the oscillation frequency ui 
over the spectrum of flux differences [5^, 



U! 



J^deuj{e) 



de 



Jo deuj{e)i^ 



(22) 



de 



de 



(16) 



The forward scattering amplitudes due to neutrinos 
and antineutrinos scattering off each other leads to a 
collective potential. 



RlJr 



2(1'..-*..) 

and B is approximately equal to t(0, 0, 1) + 0{9\^). 

The dynamics of the neutrino flavor pendulum is ap- 
proximately determined by a comparison of the potential 
and kinetic energy of the system [i^. The potential 
energy of the system is the first term on the r.h.s. of 
Eq. (|2T|) . i.e., u)T+/vr, where 



(17) 



B 



The potential weakens as at large distances because 

the fluxes dilute as 1 /r^ and there is another approxi- which is initially 



de(P 



mately 1/r suppression from the last term in brackets, 
because the neutrino flux becomes more collinear at 
large distances. The potential enters the Hamiltonian 
as Vr ^{r)T), where 



2$. 



(23) 



(24) 



D 



/ &(P(6)-P(e)) 
Jo 



(18) 



Note that the Hamiltonian now depends on P and P 
themselves, thus making the flavor evolution nonlinear. 

The single-angle equations of motion for the flavor 
composition of neutrino and antineutrino fluxes from a 
core-collapse supernova are then given by 

dP(e) / cj(e)B A(r)L , ^, , , 

^ = { + + l^{r)T> X P(e) , (19) 



is the measure of the fraction of the neutrino flux 
available for oscillation. We remind the reader that the =F 
sign depends on whether Am^j^^ > (normal hierarchy) 
or < (inverted hierarchy). On the other hand, the 
kinetic energy of the system is the second term on the 
r.h.s. of Eq. (I^TI) . given by i^(r)J^^, where 



dr 
dPje 

dr 



Vr Vr 

oj{e)B A(r)L 



which is initially 

^-(i?.J 



IDI 



(25) 



(26) 



M(r)D X P(e) . (20) 



B. Flavor Evolution 

In the central regions of a core-collapse supernova, the 
matter potential A(r) ^ w, and the mixing angle 613 is 
suppressed by the factor ^ uj/\ [H, [13] • As was shown 
in [13, m, 113], the role of a large MSW potential is 



is the net lepton asymmetry in the system. For the 
first few 100 ms after bounce, in which the explosion 
mechanism must operate, one typically has > > 
, thus P2 and P^ arc positive except at the very 
highest energies. For normal hierarchy, i.e., Ato^^jjj > 0, 
the potential energy w/w^B • de (P -I- P) is already 
negative in the initial state, and therefore the pendulum 
remains close to its initial state (any other configuration 
would have higher energy). In the inverted hierarchy, 
i.e., Am^t^j < 0, however, the potential energy is 
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positive initially, and flipping the polarization vectors 
P and P leads to a lowering of the total energy. 
This instability of the neutrino flavor distribution leads 
to almost complete flavor conversions by flipping the 
polarization vectors P and P if Aml^.^ < 0. 

The above argument assumes that all P, and Pz are 
positive, which is true for supernova neutrino fluxes in 
the early postbounce phase, at all but the very highest 
energies. This may not always be the case and it is worth 
emphasizing that even for normal hierarchy it is possible 
to get flavor flips, as first shown in [g^I- However, that 
requires a larger flux of ^ and r neutrinos (to make P^ 
and P^ negative) than is typically available in the early 
postbounce phase. 

When the kinetic energy i/i(r)J^^ becomes compara- 
ble to the potential energy ujT+/vr, the flavor pendulum 
begins flipping back and forth between the up and down 
states, and does not always return to its initial position. 
This radius, below which the collective potential pins 
all polarization vectors together and the motion is 
synchronized, is given by the condition 

M(rsync)«45(l+-^) (27) 
\ ^' sync / — 

Beyond this radius, the neutrinos begin to convert flavor 
as the flavor pendulum tends to drift towards the lower 
energy configuration. 

As the neutrinos stream out, the magnitudes 
of the vacuum Hamiltonian uj/vr and the collec- 
tive Hamiltonian fjL\D\ eventually become compara- 
ble. This happens at a radius rend where fi{rcnd) ~ 
uj [l + RlJ{4:rl^^)) the flavor pendulum settles 
into the lower energy state, which involves a flip in flavor 
space in the inverted hierarchy. Collective flavor conver- 
sions approximately freeze-out at this radius. Vacuum 
and MSW neutrino oscillations take place at much larger 
radii and we neglect them here. 

The flip of the flavor pendulum, as described above, 
leads to a swap of the i>e and i?e number fluxes with those 
of and number fluxes via pair-conversions I'ei^e ^ 
VxVx- The number fluxes after collective effects (on the 
l.h.s.) are given in terms of the number fluxes before 
collective effects (on the r.h.s.) as 

^<^^P"* ,(28) 
\d^yjde e>espiit ' 

d^pjde^ d^^Jde ,(29) 

4d$^^ /de Ide + /de + 2d$^^ /de , (30) 

where EgpUt is given by the constraint [5]| 

/ ded^i,^/de^ / de {d^i,^/de - d^pjde) . 

(31) 

The sharp discontinuity for at egput is known in 
the literature as a spectral split, and appears at 
~(6 — 8) MeV. There is a spectral split in Vf. too. 




Radius [km] 



FIG. 1: Survival probabilities of flavor states Ve (same for 
Vx) at representative energies (5, 10, 25, and 40 MeV) as 
a function of distance from the origin in a snapshot at 250 
ms after bounce in the core collapse evolution of the 15- M0 
progenitor star discussed in Section [IV Al 



but typically at even lower energies (< 5 MeV) and is 
generally ignored [H, . 

A representative example of Ve flavor evolution is 
shown in Fig.[TJ Note how the survival probability is 
initially close to one [v^ preserve their original state), 
begins to decrease from Tsync (which in the example 
shown here is ^ 250 km), and finally around Tend (here ^ 
700 km) asymptotes to zero (complete flavor conversion) . 
Note that there are neutrinos that return back to their 
original state - those are the lowest energy Vf. (below 
~ 7 MeV) that do not convert flavor as mentioned 
above. We remind that in a two-flavor approximation, 
this survival probability for is exactly the same as for 
Vx- The behavior of Vf, (and Dx) is only slightly different 
(energies below ^ 4 MeV return back to their original 
flavor), and therefore not shown. 

The value of the collective potential when 
the flavors flips start occurring, i.e., ^i{r sync) = 
4w (1 + i?^^/(4rgy„j,)) ^+/^-; docs not depend on the 
total neutrino number flux, but only on the relative 
number fluxes. Note that in ^(r) cancels with 

the denominator of and is simply a choice of 

normalization. With our normalizations for P, P, 
and /i(r), all model-dependence on neutrino spectra is 
absorbed into one number, i.e., u]F+/Ft. 

For core-collapse supernova emission parameters pre- 
dicted by typical simulations, ujJ-^/J-'^ is typically in 
the range (30 - 300) km^^ We plot in Fig.H the 
synchronization radius Tgync as a function of the number 
luminosity My^c and the neutrinosphere radius R^^ using 
Eg. (P7|) and assuming wJ"+/J"^ = 50 km~^ Note 
that a factor of six increase in the value chosen as the 
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Neutrinosphere Radius Ry^ [km] 



FIG. 2: Contours of equal synchronization radius (in km), us- 
ing Eq. (|27p . as a function of the neutrino number luminosity 
Af^,i, — i-KR^^^v^n (in 10^^ s~^) and the Ue neutrinosphere 
radius R^^ for a fiducial value of i3T+/J-l taken to be 
50 km^^. See text for more details. 



critical /i(rsync), reduces rgync by only ^ 40% because 
of the l/r^ scaling of fi. In the postbounce pre- 
explosion phase, the total neutrino number luminosity is 
AC^p ^(1 — 10) X 10^^ s~^ and i^e neutrinosphere radius 
is Ru^ ^(40 — 80) km, which leads to a typical synchro- 
nization radius rgync ^(200 — 400) km. Similarly rend is 
seen to be -(400 - 700) km. 

Our discussion so far has been based on a single-angle 
formalism. However, as we have already mentioned, 
neutrinos emitted at different angles relative to the 
neutrinosphere experience different collective and MSW 
potentials. This leads to multi-angle effects. There are 
three ways in which these multi-angle effects are known 
to become important. 

First, if the MSW potential itself is much larger 
than the collective potential weighted by the lepton 
asymmetry factor, i.e., X(r)/vr ^ /i(r)|D|, collective 
oscillations are suppressed j58| . Spelled out, this is the 
case where 

A(r) >2\/2Gf$,.p^J'-. (32) 

Second, if the v,, and Vf. fluxes are very similar, 
i.e., $1/^ « ^PeJ ^ single angle treatment is not 
appropriate due to multi-angle decoherence [s^]. In 
this case, one finds that different angular modes accrue 
random phases for both normal and inverted hierarchy. 
Thus the polarization vectors P and P shrink to zero 
due to kinematic decoherence. This can begin as soon as 
the synchronization radius is reached and leads to rapid 



flavor equilibration (all flavors assume the same spectrum 
and number flux). 

Third, if the fluxes of /i and t neutrinos and antincu- 
trinos become comparable to or greater than the or 
Ve fluxes, i.e., > ^v^, ^p^, there can be self-induced 
suppression of collective oscillations [63|. If this were to 
happen, it would delay flavor conversions to larger radii. 
However, in the flrst few 100 ms after core bounce, $1/^ is 
generally significantly smaller than <^^^ and Hence, 
self-induced suppression of collective effects, shown to 
occur when k. <^^^ k. <I>p^ [6^, is rather unlikely [ssj . 

We do not have a detailed analytical understanding 
of these multi-angle effects, and they will be studied 
numerically in the Section IIV CI Before we perform a 
more detailed numerical study, we can now provide flrst 
approximate answers to the two questions raised at the 
end of Section |TT1 

(i) What is the typical enhancement in heating that one 
can expect? It is easy to see that collective oscillations 
can give rise to almost maximal flavor conversion. All 
neutrinos and antineutrinos change their spectra, thus 
the entire Vf, and spectra can get exchanged with 
those of Vx leading to the largest possible effect that 
can be expected from any flavor changing phenomenon. 
The quantities responsible for heating, i.e., £y^(ejj and 

(e?^ ) , can get replaced by C^^ (ej^ ) , which may be 
signiflcantly higher if the luminosities in all flavors are 
comparable but energies are larger, leading to net 
enhancement of neutrino heating. 

(ii) Does this enhancement take place at sufficiently 
small radii to have a significant effect on the total 
neutrino heating! Based on typical postbounce neutrino 
emission characteristics, we expect flavor exchange to be- 
gin at rgync (200— 300) km and complete at Tend —(500— 
700) km at — lOOms after bounce. Neglectin g m ulti- 
anglc effects that could force rgync ^ 700 km [74l . [t^, 
the oscillation radii will decrease in the later postbounce 
evolution, since the neutrinospheres recede and the 
luminosities decrease with time. Hence, depending 
on the detailed dynamical evolution of a given core 
collapse event, collective oscillations may indeed occur 
at sufficiently small radii to significantly affect neutrino 
heating. In the next section, we will use full radiation- 
hydrodynamic core-collapse supernova simulations to 
obtain a more quantitative handle on the relevance of 
collective oscillations including multi-angle effects. 

IV. EFFECT OF COLLECTIVE OSCILLATIONS 
ON SUPERNOVA SHOCK REVIVAL 

To explore the potential effect of collective neu- 
trino oscillations in the core-collapse supernova envi- 
ronment more quantitatively, we perform simulations 
with VULCAN/2D, an a xisy mmetric Newtonian radiation- 
hydrodynamics code [U [73, m, [Hi. In the variant of 
VULCAN/2D that we use here, neutrino transport is han- 
dled in the multi-group flux-limited diffusion (MGFLD) 
approximation to the full Boltzmann equation, evolving 
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FIG. 3: Neutrino number luminosities (left plot) and energy luminosities (right plot) as a function of time after bounce for 
model sll.2WHW02 (left panels) and model sl5WH07 (right panels). We extract the neutrino luminosities at both the 
neutrinosphere and at 500 km. The difference in the luminosities between Rv^ and 500 km is as expected: little difference in 
the luminosities but significant differences in the and luminosities due to the large accretion luminosity, especially 
in model sl5WH07. We also show the total number luminosity, Mv,v = Mv^ + Mv^_ + , and the total energy luminosity 



the mean radiation intensity and using Bruenn's flux lim- 
iter [s^. Velocity dependence and energy-redistribution 
between neutrino groups via inelastic scattering are 
neglected. Three neutrino species are evolved, z/g and 
along with a representative /i,T flavor neutrino, v^, 
using 16 energy groups, logarithmically spaced from 2.5 
to 250 MeV. The neutrino opacities are taken from 1871 . 
More details on VULCAN/2D are provided in [13, [ll llj 

The computational grid consists of an inner quasi- 
Cartesian region that gradually transitions to an outer 
polar grid starting at a radius of 20 km and extending 
out to 5000 km with 221 logarithmically-spaced radial 
zones and 121 angular zones, covering the full 180° of the 
axisymmctric domain. Using a Cartesian center avoids 
the small time steps associated with converging angular 
zones of a polar grid. See Fig. 4 of [ssl for an example of 
our grid setup. 

It is obvious from Fig. [5] that low neutrino number 
luminosities are favorable for collective oscillations to 
occur at small radii where they may have an impact 
on neutrino heating. There is a general trend (at least 
in the ^ (10 — 2O)-M0 ZAMS mass range) for more 
massive progenitors to lead to higher postbounce neu- 
trino luminosities [l2j . Hence, in this study, we perform 
calculations with the nonrotating 11.2-M0 progenitor 
model of Woosley et al. [s^ (solar composition, model 
S11.2WHW02 in the following) and, for comparison, also 
with the nonrotating 15- A/0 progenitor model of Woosley 
& Heger [93| (also solar composition, model sl5WH07 in 
the following). 



We evolve both progenitors with the H. Shen et 
al. EOS (HShcn EOS, [H HI]), which is based on 
a relativistic mean field model of nuclear matter and 
yields a maximum gravitational mass of 2.24-M0 for a 
cold neutron star. We follow models sll.2WHW02 and 
S15WH07 from the onset of core collapse to 400 ms after 
core bounce and do not observe an onset of explosion in 
either model before we terminate our calculations. For 
the 11.2-M0 progenitor, Buras et al. [s^ observed an 
early and weak neutrino-driven, SASI-aided explosion. 
For the 15-M0 progenitor, Bruenn et al. reported an 
explosion setting in at ~ 300 ms after bounce. These 
differences in outcome observed by these groups may be 
due to their use of the softest variant of the Lattimer- 
Swesty EOS [l3| (which has now been ruled out p^). 
inclusion of general- relativistic effects and/or their more 
sophisticated treatment of neutrino transport. 

A. Postbounce Evolution: 
Neutrino Radiation Fields and Hydrodynamics 

In Fig. [21 wc show the postbounce neutrino energy lu- 
minosities {Cui ) and number luminosities (A/"^. ) for mod- 
els S11.2WHW02 and sl5WH07. We extract' the angle- 
averaged luminosities both at the neutrinosphere 
(i?j/^) and at a radius of 500 km. Significant Ve and, to 
some extent, Df, emission does occur outside the neutri- 
nosphere due to charged-current interactions involving 
accreted dissociated material [accretion luminosity] e.g., 
[93|). Since the do not participate in charged-current 
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interactions, their luminosities evolve little between R^,^ 
and 500 km. 

In model sl5WH07, the neutrino luminosities are 
consistently higher than in model sll.2WHW02. This is 
due to the higher temperatures in this progenitor, which 
lead to higher core luminosities, and to a higher accretion 
rate, which leads to higher accretion luminosities. With 
the exception of a very short period close to bounce (< 
30 ms) , the standard hierarchy of neutrino number flux 
($1,, = J\f„j4:TTr'^), $1,^ > > <I>y^, is achieved in both 
models at Ri,^ and at 500 km. Such a hierarchy is not 
present in the energy luminosities at i?^^ , but is obtained 
when taking the accretion luminosity into account. 

In the following, we make the simplifying assumption 
that neutrinos traveling through the neutrinosphere will 
undergo collective oscillations but those emitted as part 
of the accretion luminosity will not. This approximation 
is difficult to overcome, since a full collisional Boltz- 
mann solution including collective oscillations would be 
required for a self-consistent treatment. However, tests 
in which we used the asymptotic + Cf,'^'^) instead 
of the ncutrinospheric luminosities led to no qualitative 
and only small quantitative differences in the critical 
oscillation radii. 

In the left panel of Fig.|4l we show the evolution of the 
shock radii at the North and South pole, and of the angle- 
averaged shock radius in model sll.2WIIW02. Also 
shown are the energy-averaged neutrinosphere and 
gain radii in this model. Until ~^150ms after bounce, the 
shock remains essentially spherically symmetric. Then, 
the SASI begins to grow and large asymmetries arise 
in the shock front. The £ = 1 sloshing of the shock 
radius in the North-South direction is characteristic of 
the SASI in 2D (cf. M , but note that 3D gives a different 
SASI behavior, e.g., [30|)- While oscillations in the shock 
position reach radii upwards of ~400 km, the average 
peaks at ~300km at ^150 ms after bounce, then slowly 
recedes, reaching ^200 km at 350 ms. The angle-averaged 
gain radius, where charged-current neutrino heating 
balances cooling, hovers around ^100 km. The hashed 
region denotes the angle-averaged radial extent in which 
75% of the net heating occurs. The charged-current 
interactions are most effective at transferring energy to 
the matter at high density (cf. Eq. ([2])), therefore, most 
of the net heating occurs near the gain radius where 
the density is the highest (cf. the discussion in [t^ and 
Fig. 10 of [Z^l depicting the heating rate as a function of 
radius). The energy- averaged i^e neutrinosphere radius 
peaks near ~70km at ^40 ms after bounce and recedes 
thereafter. By 350 ms after bounce, the Vg neutrinosphere 
has receded to ~^40 km. 

The postbounce evolution of model sl5WII07 is sum- 
marized by the right panel of Fig.|31 The postbounce dy- 
namics is qualitatively similar to model sll.2WIIW02's 
and quantitative differences are due primarily to model 
sl5WH07's higher postbounce accretion rate. This pre- 
vents the shock from reaching the higher radii achieved 
in model sll.2WIIW02 before stagnation and suppresses 



the development and the strength of the SASI until 
later times. 



B. Collective Neutrino Oscillation Radii 

Collective neutrino oscillations may be relevant in the 
postbounce evolution if they occur within the region 
behind the shock where conditions arc conducive to 
net neutrino heating. We post-process the spectral 
neutrino fluxes predicted by the VULCAN/2D simulations 
in two ways to determine the radius at which collective 
oscillations may begin. As input to these calculations, 
we chose the neutrino spectra at the Ve neutrinosphere, 
which we average over lateral angle^. 

Our first method is using the analytic expressions of 
SectionHni We invert Eq. ^ using Eqs. (fT7| . (p2l) . (pi]) . 



and 



and solve for rgync, 



^sync Riy^ 



1 + V2Gf'^..^/{^T+/TI) 
9 + 8V2Gf $,^p/(^}^+/J-2 ) - 1 



1/2 



(33) 



We also use ^(rend) ~ uj [l + -R^^/(4r^„d)) f-^- *° obtain 
an estimate for the radius at which the oscillations 
effectively are complete. Once again inverting this using 
Eqs. (Ull), (HH), and we obtain for rend, 



^cnd — 



1 + AV2Gf'S?^.pJ'-/^ 
9 + 32V2GF'i>u,uJ'-/^ -I 



1/2 



(34) 



As an alternative to the rather rough approximation 
of Eqs. (p3| and (p4)) . we determine the critical collective 
neutrino oscillation radii by numerically solving the set 
of coupled nonlinear differential Eqs. ([T9l) and ([20|) with 
the initial conditions given by Eqs. ([S]) and ^ based 
on the neutrino number fluxes at R,y^ in the VULCAN/2D 
simulations. The equations are solved as a function of 
radius r for 32, 64, and 128 energy groups spaced as in 
Gauss-Lcgendre quadrature with the oscillation code of 
Dasgupta et al. 5J] in combination with the open-source 
ordinary differential equation solver CVODE [95|. We 
carry out this calculation for select postbounce times 
and numerically identify rgync and rend with the radii 
at which 5% and 90% flavor conversion have occurred, 
respectively. For reference, Fig.[T] shows the survival 
probability for select energies obtained with such an 
evolution for the neutrino spectra in model sl5WH07 at 
250 ms after bounce. 

We present the results of both methods applied to 
models sll.2WHW02 and sl5WH07 in Fig.H The 



values of 



sync 



and Tend predicted by the two methods 



agree well at early to intermediate times. At late 



^ Since our models are nonrotating, the variation of the neutrino 
spectra with lateral angle is not large (cf. [M [Z3i [Hi H ) • 
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FIG. 4: Comparison of neutrino oscillation radii to shock and gain radii. Average (thin solid lines) and polar shock radii 
(dot-dot-dashed and dot-dashed-dashed lines) in the postbounce phase of model sll.2WHW02 (left plot) and model sl5WH07 
(right plot). Also shown (dot-dashed lines) are the energy-averaged i^e neutrinosphere locations and the angle- averaged gain 
radii (dashed lines). The hashed region just above the gain radius denotes the area in which 75% of the net heating occurs. We 
show the radii at which collective neutrino oscillations begin (rsync; dark black lines and symbols) and end (rond; light green 
lines and symbols) as obtained via an analytic approximation (solid lines) and via detailed numerical calculations at select 
times (squares and circles), rsync is initially large, but drops to radii comparable to the shock radii at ~(100— 150) ms in model 
S11.2WHW02 and at ~(300 - 350) ms in model sl5WH07. See text for details. 



times, the analytic approximation ovcrpredicts by ^25%. 
Note that, although the analytical formulae derived in 
Sec. lIIIBl are in a two-flavor approximation, they can be 
compared to a three-flavor numerical calculation, because 
the third flavor is almost decoupled. 

It is apparent from Fig. 2] that the radial interval over 
which the collective oscillations occur is well outside the 
shock at early times. Only after ^^ISOms (~350ms) 
in model sll.2WHW02 (sl5WH07), when the total 
number luminosity has decreased and the neutrinosphere 
has receded, does rgync recede below the shock radius 
and collective oscillations can have an effect on the 
subsequent evolution. 

C. Multi-angle Effects 

To ascertain the importance of multi-angle effects 
we use the neutrino luminosities and electron density 
profile as predicted along the equatorial direction by our 
VULCAN/2D simulations. Due to the high computational 
demand of multi-angle oscillation calculations, we do 
not include the energy spectra of the different flavors, 
and instead assume a monoenergetic ensemble with 
the average vacuum oscillation frequency w. This is 
the same approximation as used by [TJ, [z3- Since 
our VULCAN/2D simulations made use of the efficient 
MGFLD variant of the code, which does not carry direct 
information on the momentum-space angle dependence 



of the neutrino radiation field, we compare with the 
fully angle-dependent calculations of [t^I and construct 
approximate angle- dependent radiation fields. We find 
that the angular distribution of the neutrino luminosity, 



derived from the simulations of 77[, is parametrized 
quite well by d^/dcosdji^^ oc exp [(cosi^^^^ — l)/a'_R„J , 
where we choose a^^ = 0.357((1/J-") — 1), and where 
i^/J^) is the inverse flux factor. This parameterization 
reproduces both the isotropic ^ 1, a ^ 1, and 

d^/dcos-d oc constant) and the free streaming ~ 
1, cr '--^ 0, and d^/dcosi? cx 5{'&)) limits, and qualitatively 
reproduces the angular distribution at the neutrinosphere 
(cf. Figure 3 of |77l|). At the neutrinosphere, our inverse 
flux factors are ^4-5. We assume a sharp neutrinosphere, 
cutting off all neutrinos traveling backwards into the 
neutrinosphere. With these choices, we solve Eqs. ([13]) 
and (HH) using a multi-angle oscillation code that is 
technically similar to the one used for the single-angle 
calculations.^. 



^ We have verified numerical convergence by comparing results 
from calculations with 100, 200, and 400 angular bins. The local 
error tolerance was fixed at 10~^^, which allows us to achieve 
convergence with ~ 400 modes. Additionally, we have verified 
the calculations with 800 modes in a limited number of cases, and 
found them to be consistent. We have also reproduced results 
similar to Ref. l52l . 
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FIG. 5: Survival probabilities (averaged over all emission angles) calculated using a multi-angle oscillation code for a single 
average energy. The different panels correspond to different times for both models - t = 150 ms (left), t — 250ms (center), and 
t = 350 ms (right). Blue solid lines denote results for model sl5WH07 and dashed red lines are results for model sll.2WHW02. 
We also plot results of our cross-checks: dash-dotted green lines are results for model sll.2WHW02 with half-isotropic angular 
emission and dot-dot-dashed black lines show for the same model the case with half-isotropic angular emission and 10 times 
higher electron density. The vertical lines are from our single angle calculations and represent the radii at which 5% flavor 
conversion has occurred in models sll.2WHW02 (dotted) and sl5WH07 (dashed). 



In Fig. [SJ wc show the survival probability (averaged 
over all emission angles) of the average-energy neutrinos 
calculated using the multi-angle code. The three panels 
show the situation at 150ms (left), 250ms (center), and 
350 ms (right) after core bounce and results for model 
S11.2WHW02 and sl5WH07 are shown in dashed red 
and solid blue lines, respectively. The vertical lines in 
Fig.[5]are from our single angle calculations and represent 
the radii at which 5% flavor conversion has occurred for 
the S11.2WHW02 (dotted) and the sl5WH07 (dashed) 
progenitors. 

We find that in model sll.2WHW02 the onset of 
flavor conversion is delayed by multi-angle effects to 
~ (300 — 500) km. This is expected and due to high 
electron number density. As we remarked in Section lllll 
whenever the MSW potential A(r) ^ Ama(^) = 
2^/2Gi.$^,p(i?2^/r2)J■_, muh i-anglc effects arc strong 
and suppress the oscillations. In Fig. [51 we plot A(r) 
and Ama(?') as a function of radius at 150 ms (left 
panels), 250ms (center panels), and 350ms (right panels) 
after bounce, for models sll.2WHW02 (top panels) and 
S15WH07 (bottom panels). Radial profiles along ten 
lateral directions from the North to South pole are shown 
to capture variations in A(r) due to SASI oscillations. 
A(r) is almost always larger than Ama, which decreases 
(X as expected. A(r), which effectively traces the 
electron number density, falls off cx below the shock 
radius. Above the shock, where matter is essentially in 
free fall, the density roughly follows r~^-^. Typically, 
the ratio X{r) / XuAir) is in the range (1 — 100), getting 
close to 1 at later times (i > 250 ms) at r ~200km in 
model S11.2WHW02. 

Comparing A(r) along different directions, we find 
that the SASI oscillations developing at t — tbouncc ^ 
150 ms lead to a significant spread in the values of A(r), 



sometimes, in model sll.2WHW02, bringing it just below 
the critical value below which multi-angle suppression 
is less effective. This occurs when the shock recedes 
to particularly small radii and the region in which the 
suppression is lifted is always at or outside the shock. 

On the other hand, as we remarked in Section IIIIBl 
multi-angle effects may appear in a second way: If the 
i^e and De fluxes arc very similar, i.e., <f>^^ ~ 
multi-angle decoherence of flavors sets in. In our model 
S11.2WIIW02. the i^e and Dg fluxes are not too similar 
((<i>^^ — ^i^J/^p^ w 0.25), thus the decoherence effect 
remains negligible. 

In model sl5WII07, on the other hand, the multi- 
angle oscillation calculation does not predict a signiflcant 
delay in radius compared to the single-angle prediction, 
in spite of the high electron number density also present 
in this model. Since in model sl5WH07 the ly^ and Dg 
fluxes are quite similar (($i/^ — <f>pj/$p^ w 0.1) early on, 
we attribute this surprising result to (at least partial) 
multi-angle flavor decoherence, which leads to oscillations 
despite the high electron number density. We flnd flavor 
conversion almost as soon as the neutrinos cross the 
synchronization radius. This observation is consistent 
with the previous work by Esteban-Pretel et ai, in which 
where they found that {^^^ — ^ 0.2 leads to 

decoherence fs^ . 

Comparing our results quantitatively with the recent 
multi-angle work of Chakraborty ct al. [zl, [tH is not 
straightforward, since these authors based their calcu- 
lations on ID supernova simulations using a different 
progenitor model (a 10.8-Mq progenitor of [s^). How- 
ever, we note that they observe a stronger multi-angle 
suppression than borne out by our models. For example, 
flavor conversion is delayed for their 10.8-Mq progenitor 
almost until (700 — lOOO) km, while we observe flavor 
evolution to begin around already at (300 — 500) km in 
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FIG. 6: The MSW potential A(r) along various directions (thin solid lines, 10 rays equally spaced in cos(6'iat)), in comparison to 
the minimum A(r) needed for multi-angle suppression Ama = 2^2GF^v,i?{Fli^/r^)J'- (thick dot-dashed line, taken along the 
North pole (NP), South pole (SP), and equator (EQ)). To guide the eye, rulers show r~^'^, r~^, and r~'^ radial profiles. The 
steep rise in the A(r) profiles occurring around r ~200 km, but varying with lateral angle, is the location of the shock. The MSW 
potential A(r) is generally larger than the critical value AMA(r) needed for multi-angle suppression of neutrino oscillations. 



our models. We believe that this is due to two reasons: 

Firstly, a half-isotropic angular distribution 
(d^Jd cos oc cos ) was used in the calculations 
of [7^ [rsj . Compared to our angular distribution, 
this underestimates neutrinos emitted at large angles. 
The collective interaction is stronger for tangentially 
emitted neutrinos than for radially emitted neutrinos, 
therefore suppressing the tangential modes leads to a 
slower growth of the collective instability. We verify this 
claim by replacing the angular distribution in our model 
sll.2WHW02 by a half-isotropic angular distribution. 
The results for this half-isotropic case are shown in 
Fig. [5] (dot-dashed green line). The "double-step" 
feature only appears when we use our angular emission 
spectrum, which introduces stronger collective effects. 
This suggests that the feature is related to the angular 
spectrum, and not a numerical artifact. It is thus 
clear that simply changing the angular distribution 
can change the onset of flavor conversion. This may 
be strong enough to create a qualitative difference as 
evident from the snapshot sX t = 350 ms in our 11.2-A'f0 
model, where we find that the oscillations do not occur 
for a half-isotropic angular distribution, but do occur 
for our angular distribution modeled after the full 2D 
multi-angle neutrino transport simulations of (77j . 

Secondly, in the 10. 8- M© model of [zl, [zl, one 



finds a ratio of electron to neutrino density that is 
up to 10 times larger than in our 11.2-Mq model at 
various radii and times. This is due to the different 
progenitor structure used - our sll.2WHW02 model has 
a lower postbounce accretion rate than the models of 
[tS [75| . We verify that this is indeed an important 
factor, by artificially increasing the electron density by 
a factor of 10 and replacing the angular distribution in 
our oscillation calculations by a half-isotropic angular 
spectrum for model sll.2WHW02. This case is expected 
to closely follow the results of [tI, [t^ . The results of 
this are shown in Fig. [5] (dot-dot-dashed black line). 
They demonstrate that such a change in the electron 
density and angular distribution can indeed significantly 
suppress the flavor evolution, in agreement with previous 
results [tIIiI- 

In the light of these results, we believe that the role 
of multi-angle effects remains an issue that requires 
further scrutiny. The role of the matter density, Ve/i^e 
asymmetry, and the angular distribution need to be 
studied in more detail. Flavor conversion are not 
always completely suppressed due to multi-angle matter 
suppression. Predictions of the neutrino flavor content 
at early times must therefore be used with abundant 
caution. Fortunately, as we shall see in the next 
section, our conclusions regarding the impact of collective 
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oscillations on the supernova mechanism remain largely 
unchanged. 



D. Potential Enhancement of Neutrino Heating 



As discussed in Section Hill collective neutrino oscilla- 
tions (in the inverted mass hierarchy) will lead to a swap 
of the i^e and spectral fluxes with the spectral fluxes of 
the neutrinos (cf. Eqs. In order to illustrate 

the potential enhancement of neutrino heating due to this 
swap, we first consider the original heating rate before 
oscillations (cf. Eq. ([3])), which can be expressed as 



CP ; 



(35) 



where we explicitly split the luminosity into core lu- 
minosity emanating from the neutrinosphere (>C"^) and 
accretion luminosity {Cf,'^'^), the latter being emitted 
almost entirely interior to the gain region. We take (e^.) 
as the value calculated at the Vg neutrinosphere. The 
do not take part in charge-current interactions and play 
no role in the heating before oscillations. 

We now estimate the heating rate after taking into 
account a partial conversion of the neutrino spectra in 
the region behind the shock, as observed at late times in 
our simulations. 
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The stars {*) on (e^ ) and {e^^ ) denote that the respective 
averages are taken over only the part of the spectrum 
below the split energy egpiit (cf. Section Hill . and Cp 
are the column number densities of interacting baryons 
taking oscillations into account, 



dr Hi Pc^{r) , 



(37) 



where -Pcx(^) is the flavor conversion fraction and Rg and 
Rs are the gain and shock radius, respectively. If no 
oscillations occur, Pgx = everywhere and cf is zero, 
Rafter = ^before- If Complete oscillations occur before 
Rg7 ^ox = 1 everywhere and cf = Ci. All other quantities 
in Eq. ([55)) have their pre-oscillated values. 

We employ the neutrino data from models 
S11.2WHW02 and sl5WH07 and the analytic 
approximations to the heating, Eqs. ([55]) and ([M]). 
to calculate the change in the heating rate due to 
collective oscillations. The column number densities, Ci 
and cp, going into the heating rates are angle-averaged 
values obtained from simulation data. The results are 
depicted by Fig.[71 Shown with dashed lines are the 



1000 



^ 100 



a 

o 



10 



O.I- 



rs 

X 0.01 



0.001 



si 1.2WHW02 - expected (single-angle) 
S11.2WHW02 - expected (multi-angle) 
s 1 1 .2WHW02 - r,„j < Tj,,,,, (hypothetical) 
S15WH07 - expected (single-angle) 
S15WH07 - expected (multi-angle) 
S15WH07 - rj„j < r^,,, (hypothetical) 




0.15 0.2 0.25 

t ~ t-bounce [s] 



0.4 



FIG. 7: Time evolution of the potential percentage increase 
in the heating rate, 'Hafter/'Hbcforc — 1, due to collective 
neutrino oscillation. The dashed lines (thick for model 
S11.2WHW02 and thin for model sl5WH07) assume the 
hypothetical case of complete conversion already below the 
gain radius, leading to an enhancement of (20 — 100)% 
depending on the progenitor and time after bounce. In 
our simulations, complete conversion does not occur before 
the gain radius. The more realistic estimate of the heating 
enhancement based on the oscillation calculations in Section 
mil is much lower and shown in solid lines. Before t = 150 ms 
and 330 ms the synchronization radius is outside the shock in 
model S11.2WHW02 and sl5WH07, respectively. The points, 
blue squares for the sll.2WHW02 model and red circles for 
the sl5WH07 model, represent our estimate of the heating 
enhancement if the multi-angle survival probabilities are used. 
At 150 ms for the sll.2WHW02 model and at 150 and 250 ms 
for the sl5WH07 model, no heating enhancement is seen, the 
conversion occurs completely outside the shock. 



expected heating enhancements in model sll.2WIIW02 
(thick lines) and model sl5WH07 (thin lines) in a 
hypothetical scenario in which we assume complete flavor 
conversion below the gain radius (Pox = 1 everywhere). 
In this extreme case, the heating would be enhanced 
by >60% in model sll.2WHW02 and >20% in model 
S15WH07, this is similar to the configuration of Suwa et 
al. (TOj . Note, however, that our assumption that 
neutrinos of the accretion luminosity do not undergo 
oscillations may be invalid in this hypothetical situation. 

For a more realistic estimate of the heating enhance- 
ment, we use Tsync as numerically computed for both 
models in Section llV Bl and assume that above this radius 



0.05 and below this radius, P,- 



0. As a 



consequence, if rgync is greater than the shock radius, 
there is no heating enhancement. The result of this 
is shown in Fig. [7] in thick-solid (thin-solid) lines for 
model S11.2WHW02 (sl5WH07). The predicted heating 
enhancement sets in at much later times and is < 0.1% for 
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both models. If one instead assumed that Pex = 1 outside 
of ''sync, the enhancement would still be < (2 — 3)% 

When considering the predicted flavor conversion in 
our multi-angle calculations discussed in Section llV C| we 
find, as expected, that there is no further enhancement 
of the heating, but rather that the enhancement is 
even more suppressed. We show this in Fig. [71 where 
we denote by points the expected heating enhancement 
using the multi-angle survival probabilities presented in 
Fig. |6] at three postbounce times for each progenitor. 
In the S11.2WHW02 model, the heating enhancement is 
zero or reduced significantly. The heating enhancement 
estimated in the sl5WH07 model, where we observed 
little change in the onset of collective neutrino oscilla- 
tions, is not as strongly affected. Hence, wc conclude 
that collective neutrino oscillations are very unlikely to 
have a measurable effect on the neutrino heating and 
postbounce dynamics in progenitors in and above the 
mass range considered in this study. 

V. DISCUSSION AND CONCLUSIONS 

Almost eight decades after Baade & Zwicky's stageset- 
ting 1934 proposal [9^,|93] that a core-collapse supernova 
represents the transition of an ordinary star to a neutron 
star, the details of the mediating mechanism that 
converts gravitational energy of collapse into energy of 
the core-collapse supernova explosion remain uncertain. 

The neutrino mechanism, based on net heating by 
charged-current neutrino absorption in the region just 
below the stalled shock, appears to be the most vi- 
able candidate mechanism, requiring the least spe- 
cial conditions {e.g., not requiring rapid rotation or 
strong magnetic fields etc.) to succeed in exploding 
garden-variety Type-II supernova progenitor stars. Yet, 
neutrino-driven explosions fail in ID, are marginal in 
2D simulations, and modeling groups are now exploring 
the neutrino mechanism's potentially boosted efficacy 
in 3D [2^ [s^]. While dimensionality may be the key 
to successful explosions, it is also possible that current 
1D/2D models are still missing some physics crucial to 
successful explosions. 

In this paper, we have considered new physics pre- 
viously left out of core-collapse supernova models: col- 
lective neutrino flavor oscillations induced by neutrino- 
neutrino forward scattering in the core-collapse super- 
nova core. If the neutrino mass hierarchy is inverted 
(Am^jj,^ < 0), collective oscillations will invariably 
lead to a swap of Ve and Pe spectra with the signif- 
icantly harder spectra of their heavy-lepton neutrino 
counterparts. Assuming a hypothetical scenario in which 
this swap is complete below the gain region, we find 
that neutrino heating is enhanced by (20 — 100)% in 
representative Type-II supernova progenitors of 11.2-Mq 
and 15-Mq. Such a significant boost of heating may lead 
to strong, early explosions, breaking the strong feedback 
between EOS, weak interactions, neutrino transport, and 
hydrodynamics that is present in the postbounce phase of 



core-collapse supernovae and that tends to absorb small 
variations in any of its components'^ . 

To study the viability of this scenario, we have per- 
formed collective neutrino oscillation calculations in the 
single-angle and multi-angle approximation on the basis 
of the neutrino radiation fields obtained from 2D neutrino 
radiation-hydrodynamic core-collapsc simulations using 
II.2-M0 and 15-M0 progenitors. From the oscillation 
calculations we obtain the characteristic radii 7'sync and 
Tend at which ^5% and ^--^90% of the flavor conversion 
have occurred, respectively. In our simulations, these 
radii start hundreds of kilometers above the typical shock 
radii in the early postbounce phase, but recede with 
time as the neutrinosphcres settle to smaller radii as the 
neutrino luminosities decrease. The radius of onset of 
oscillations (rgync) reaches the average shock radius at 
(150 — 350) ms and thereafter stays close to the latter, 
while Tend is systematically (150 — 200) km outside the 
average shock radius. As a consequence, most of the 
flavor conversion occurs outside the shock and can have 
no effect on the postbounce heating and hydrodynamics. 
Those oscillations that take place inside the shock occur 
close to the shock radius and, even when taking large 
shock excursion driven by the SASI into account, enhance 
the net heating by less than ^^(2 — 3)% in the most 
optimistic case. These results strongly suggest that 
collective neutrino oscillations are unlikely to have any 
qualitative or significant quantitative impact on the 
postbounce evolution and the explosion mechanism in the 
standard Type-II supernova progenitors considered here. 

Our results also show that the characteristic oscillation 
radii assume small values faster in cooler, less massive 
progenitors with lower postbounce accretion rates and 
smaller neutrinosphere radii in the postbounce phase. 
While oscillations are unlikely to boost the pre-explosion 
neutrino heating in our 11.2-AfQ and IB-Mq models, the 
situation may be different in even lower-mass progenitors 
with 0-Ne cores or in 0-Ne white-dwarf progenitors of 
accretion induced collapse. The weak explosions already 
obtained for such progenitors [l-H, 19^ could thus be 
significantly enhanced. 

As an aside, we have also shown that the multi-angle 
matter suppression [H, [zl, [t^ is somewhat sensitive to 
the choice of angular emission spectrum and the matter 
density profile. This is expected to affect the predictions 
of flavor evolution. 

In deriving our results, we have made some ap- 
proximations and simplifications. Of these the most 
limiting is that we have assumed a sharp neutrinosphere, 
common to all flavors and energies. This ignores the 
interplay of oscillations and collisions, but a completely 
self-consistent treatment would require a full collisional 
Boltzmann calculation including oscillations, which must 
be left to future work. Also, due to computational 



^ This behavior of strongly coupled astrophysical systems is known 
as Mazurek's Law |98|| . 
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limitations, we have not carried out multi-energy multi- 
angle oscillation calculations, which could be improved 
upon in subsequent work. The other significant limita- 
tion is the assumption of axisymmetry in our supernova 
calculations. Future, full 3D radiation-hydrodynamics 
simulations may lead to different hydrodynamic post- 
bounce evolutions and could yield different results. Given 
that the understanding of core-collapse supernova physics 
and collective neutrino oscillations is still in a state 
of rapid development, there might be additional, yet 
unknown, effects that could change our conclusions. 

Finally, to summarize, collective neutrino oscillations 
are an intriguing phenomenon. They have important 
important ramifications for the core-collapse supernova 
neutrino signature, but, as we have shown in this paper, 
they develop at too large radii to play a significant role 
in the explosion mechanism. 
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